Tr. J. of Physics 
XX (2002) ,1 - 8. 
© TUBITAK 



On complex permittivity of dilute random binary 
dielectric mixtures in two-dimensions 

Enis Tuncer * 

Chalmers University of Technology, SE-^12 96 Gothenburg Sweden 

Emre Tuncer 

Monterey Design, 894 Ross Drive, Sunnyvale CA 94089 USA 
Received February 1, 2008 



Abstract 

Influence of number of particles considered in numerical simulations on complex 
dielectric permittivity of binary dilute dielectric mixtures in two-dimensions are 
reported. In the simulations, dodecagons (polygons with 12-sides) were used to 
mimic disk-shaped inclusions. Using such an approach we were able to consider 
16 particles in a unit-square. The effective dielectric permittivity of the mixtures 
were calculated using the finite element at two different frequecies which were much 
higher and lower than the characteristic relaxation rate of the Maxwell- Wagner- 
Sillars polarization. The results were compared to an analytical solution. It was 
found that the permittivity values at low frequencies were inside Wiener bounds 
however they violated the Hashin-Shtrikman bounds. 

Key words: Dielectric permittivity, the finite element method, effective medium 
theory. 

The complex permittivity of materials needs to be estimated in many applications [1, 
2, 3, 4, 5]. By using classical methods, i.e., analytical techniques, effective medium 
approximation, etc., the determination of the complex permittivity of dielectric materials 
is usually carried out for simple structures [1]. Moreover, the computed values should be 
verified by experiments or by means of other computational methods. Recently, the finite 
element (fe) computer simulations [6] has shown that disk-shaped inclusions in binary 
mixtures can be imitated by n-gons (polygon with n-sides) such that when n > 10 the 
obtained electrical properties of the mixture were within 0.1 % error of the analytical 
solutions when the mixture was a dilute solution. In this paper, we investigate the 
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electrical properties of a dilute binary mixture by considering the number of particles in 
the computational domain and compare the results with those obtained analytically [7]. 

In order to obtain random structures, an algorithm based on acception and rejection 
method (also known as random sequential adsorption) was constructed [8]. In this ap- 
proach, hard-disks with radii, r, were generated with random center coordinates in a unit 
square. Overlapping of disks both with outer boundaries of the unit square and with 
other disks were not allowed. The coordinate-values x and y were picked from a uniform 
random number generator where (r, r) < (x, y) < (1 — r, 1 — r). The first hard-disk was 
fixed, then, the number of hard-disks was incrimented. In each step, the condition for 
overlapping of the new generated disk with the already fixed disk(s) was checked. In 
case of overlapping, the disk was rejected, and a new one was picked from a new se- 
ries of coordinates. In the process, obviously both the concentration of the hard-disks 
in the background medium and the number of disks should be kept constant. As men- 
tioned previously, instead of using exact disk shapes and consuming lot of computational 
power, dodecagon (polygon with 12 sides) were used in the simulations. The considered 
dodecagons were monodispersed. The boundary conditions in the fe were such that 
reflections of the computational domain in the 4 facing-directions, von Neuman neigh- 
bors [9] were used which was also Rayleigh's approach [10]. The fe simulation procedure 
has been presented previously elsewhere [11, 12, 6]. 

In the simulations the concentration of the inclusion phase q and number of inclusions 
in the computional domain were considered as variables. Only four concentration levels 
were included, q = {0.1, 0.2, 0.3, 0.4}. Starting with one inclusion, we have increased 
the number of inclusions n, in the computations. The meshing of the region was arranged 
by the sizes of the dodecagons such that smallest meshing element was 1/4 of the length 
of the dodecagon side. This kind of procedure has allowed us to have enough number of 
discritization points in the inclusions. Computations with larger number of inclusions has 
increased the number of mesh elements, therefore, we were not able to obtain solutions 
for n > 16 2 and q > 0.2. For the q — 0.1 case, the number of mesh elements were high 
enough when n > 12 2 . 

The electrical properties of the constituent media without any intrinsic polarizations 
were choosen such that the relaxation rate of the interfacial polarization in the composite 
medium [13, 14, 15] was around 1 Hz. The complex dielectric permittivity of a medium 
as a function of angular frequency, u) — can be expressed as, 

£i{uj) = e t + (1) 
ie ui 

where Si is the dielectric permittivity of the medium at optical frequencies and cr^ is the 
conductivity of the constituent, eo is the permittivity of the free space, e = 1/36-7T nF/m 
and i = \f— T. The permittivity and conductivity of the matrix phase, e\ and <j\, were 
taken 2 and 1 pS/m, respectively. The permittivity and conductivity of the dodecagon in- 
clusions, e 2 and (72, were 10 and 100 pS/m, respectively. We have estimated the electrical 
properties em, £^ and a where subscripts M and \x represents the frequency. 

In Fig. 1, the electric potential and electric field distributions of one of the considered 
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(a) (b) 

Figure 1. Potential (horizantal lines) and electric field (vertical lines) distributions at (a) lower 
and (b) higher frequencies than the characteristic interfacial polarization is observed. Number 
of inclusions are 36 and the concentration q is 0.2. 



case with n — 6 2 and q = 0.2 are presented. The difference between the two pictures are 
the one on the left is at low frequencies [y = 1 uHz) where the field distribution in the 
material is determined by the resistive properties of the media. Similarly, the picture on 
the right is at v = 1 MHz where the capacitive properties of the media set the effective 
properties of the mixture. 

The Maxwell- Garnett dielectric mixture formula [16] is valid for two-dimensional di- 
lute mixtures when the shape of the inclusions are disks [6]. The effective complex 
dielectric permittivity of the mixture, e e {to) can be found from 



e e (ui) + (n - l)ei(w) £2^) + (n - l)ei(w) 

where n = 2 in two-dimensions [12]. In Fig. 2, calculated electrical properties are presented 
as normalized to the values obtained from Eq. (2). In the figure, the bounds for dielectric 
permittivity and conductivity obtained from Wiener [17] and Hashin-Shtrikman [18] are 

also displayed as dashed (- - -) and solid ( ) lines respectively. The Wiener bounds 

are obtained when n = and n = 00 in Eq. (2). In two-dimensions, the Hashin-Shtrikman 
bounds are expressed as 

~ q 1 - q . . 

£l (£2 - si)- 1 + (1 - 9)(2ei)- 1 - £ - 62 + (ei - S2)- 1 + g(2e 2 )- 1 { ' 

The conductivity of the mixture was calculated from the dielectric losses, a — 3[e" (uj)eouj] 
at lo = 2ir uHz. It is significant from Fig. 2 that the lower bound of Hashin-Shtrikman 
is for disks since the obtained values ratios of the bounds and Eq. (2) with n = 2 are 
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Figure 2. Ratio of electrical properties calculated and estimates of the analytical formula as a 
function of inclusion concentration q; (a) e M , (b) e H and (c) a. The dashed (- - -) and solid 
( ) lines represent the bounds of Wiener and Hashin-Shtrikman. 
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Figure 3. Ratio of electrical properties calculated and estimates of the analytical formula as 
a function of square root number of inclusions n; (a) em, (b) e H and (c) a. Different symbols 
represent the considered inclusion concentrations q; (Q) 0.1, (□) 0.2, (O) 0.3 and (A) 0.4. 
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equal. This proves that the considered dilute binary mixture was isotropic. The error 
bars in the figure or in other words the spread of the electrical properties indicate the 
importance of local field effects and interaction between inclusion particles such that at 
high concentrations the interaction between particles started to influence the electrical 
properties. 

When e at low frequencies are taken into consideration (Fig. 2a), the upper limits of 
the two bounds have shown extraordinary behavior for higher concentrations of inclusions, 
the validity of Eq. (3) was very narrow, < q < 0.2. The Wiener bounds are on the 
other hand valid for our simulations, however, for q > ir/4 the Wiener bounds also show 
extraordinary behavior. It is remarkable that q — 7r/4 is the limiting concentration for 
regular arrangement of inclusions on a square net. 

The permittivity values at high frequencies are displayed in Fig. 2b, and they were 
similar to the analytical formula in Eq. (2), and the spread of the data was very small 
compared to the conductivity values presented in Fig. 2c. This was due to the ratio 
between the permittivity values considered, 3?[ei]/3?[e2] < 3[ei]/3[?2] [19]. For high 
frequency permittivity and ohmic conductivity both bounds were valid. 

When the influence of the number of inclusions in the computational domain arc 
considered, it was observed that the spread of the electrical values decreased as the 
number of inclusions increased. This behavior is displayed in Fig. 3 again as the ratio 
between the computed values and Eq. (2). It was clear that the error-bars decreased by 
increasing the number of inclusions. There was a jump between the electrical properties 
when one-inclusion and four-inclusions systems were compared to each other. This was 
due to the change in the field distributions around inclusions which has pointed out 
the importance of boundary conditions in the simulations. At small scale systems with 
less number of inclusions, the boundary conditions are crucial since each particle is not 
affected by the local field distributions created by the neighboring inclusions. In addition, 
by increasing the number of inclusions in the computation domain, the probability of 
having larger number of neighbors in a radius is also increased when compared with a 
regular arrangement of inclusions [12]. This increase in the number of neighbors should 
in principle enhance the polarization of individual inclusions which is more realistic to 
real materials. This approach was desired but not selected due to difficulties arising in 
the field calculations, mainly the meshing procedure in the fe method. The distances 
to nearest neighbors were not influenced by the number of inclusions if scaled with the 
inclusion radius. One could obviously use as many hard-disks as possible, after all it 
would be the most realistic condition. This could be done as proposed by considering 
a dodecagon instead of a circle. Finally, similar approach sould be applied to three 
dimensional problems to increase the number of considered particles in simulations. 
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